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1.0  INTRODUCTION 

High  power  microwave  amplifier  tubes,  Including  Injected  beam  crossed 
field  amplifiers  (IBCFA's)  and  broad  band  traveling  wave  tubes  (TWT's)  may 
exhibit  several  types  of  spurious  signals  In  their  output  power  spectrum. 

This  paper  describes  one  class  of  spurious  signals,  called  parametric  oscil¬ 
lations,  for  which  a  mathematical  model  Is  proposed. 

By  charting  the  frequency  f$p  of  the  parametric  signals  as  a  function 
of  drive  frequency  f^,  an  empirical  relationship  Is  obtained: 

fsp  *  mf o  -  nfd  0) 

where  n  and  m  are  positive  or  negative  Integers.  f0  Is  a  fixed  frequency 
which  Is  correlated  empirically  to  a  limit  frequency  In  the  stopband  of  the 
slow-wave  structure.  Figures  1  and  2  show  typical  measured  data  obtained 
from  an  S-band  IBCFA  and  a  Ku-band  THT. 

Relation  (1)  seems  to  Indicate  that  the  origin  of  these  oscillations 
Is  due  to  Intermodulatlon  products  between  two  signals  In  the  amplifier. 

Assuming  that  the  spurious  signals  are  Intermodulatlon  products  between  the 
signal  frequency  and  an  oscillation  Inherent  In  the  tube  at  the  frequency  f0, 
then  one  would  expect  to  find  a  coherent  oscillation  In  the  absence  of  the 
Injected  signal  (drive).  Such  Is  not  the  case.  In  fact,  high  drive  power 
levels  are  required  for  these  oscillations  to  occur.  Figure  3  Illustrates 
this  behavior.  Here,  the  amplifier  Is  driven  at  a  frequency  of  3.9GHz,  and  the 
lower  edge  of  the  circuit  stopband  Is  known  to  occur  at  7.8GHz.  The  upper 
right  picture  of  the  spectrum  analyzer  display  (Figure  3)  shows  a  normal  clean 
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output  power  spectrum  from  the  amplifier  driven  close  to  (or  just  Into) 
saturation  (drive  power  level  32W).  The  output  power  Is  measured  to  be 
1 . 6k  W. 

The  photograph  In  the  upper  left  hand  corner  of  Figure  3  shows  the 
spectrum  containing  the  spurious  signals  when  the  amplifier  Is  driven  Into 
hard  saturation  (drive  power  level  s  45W).  The  photograph  In  the  lower 
left  comer  of  Figure  3  shows  the  output  spectrum  appearing  under  these  same 
conditions  at  the  harmonic  frequency  of  7.8GHz.  The  output  power  dropped 
from  1.6  to  lkH  at  the  signal  frequency,  resulting  In  a  " power- suckout" 
frequently  observed  In  TWT's  and  CFAls. 

Since  the  spurious  signals  are  only  observed  under  full  beam  power  con¬ 
ditions  and  with  the  presence  of  a  drive  signal,  preferably  at  high  power. 

It  Is  believed  that  the  amplifier  Is  pumped  by  a  signal  related  to  the  drive 
signal,  similar  to  a  parametric  amplifier.  Hence,  we  have  coined  the  term 
"parametric  oscillations"  to  describe  these  spurious  signals. 

A  model  which  describes  parametric  oscillations  must  be  based  on  the 
existence  of  a)  a  full  power  electron  beam,  b)  a  strong  RF  drive  signal,  and 
c)  a  periodic  structure. 

A  model,  fulfilling  these  three  requirements,  was  proposed  to  the  Air 
Force  Office  of  Scientific  Research  for  further  evaluation,  and  sponsored 
by  the  agency  from  July  1977  to  June  1979.  The  first  phase  of  the  study 
program  applied  the  model  to  Injected  beam  crossed  field  amplifiers.  The 
model  predicted  the  simultaneous  existence  of  very  low  frequency  oscillations, 
which  have  been  detected  experimentally.  These  results,  which  have  been 
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discussed  In  detail  In  an  Interim  report  submitted  to  AFOSR^),  In  pub¬ 
lication2)  and  at  the  IEDM  meeting  In  Washington,  19783),  will  be  summarized 
In  this  report. 

The  second  phase  of  the  program  applied  the  model  to  helix  type  traveling 
wave  tubes.  The  results  obtained  will  also  be  presented  in  this  final  report. 


^  6.  Dohler,  Parametric  Oscillations  In  High  Power  Microwave  Amplifiers, 
Contract  No.  F49620-77-C-0096  (1979). 

2)  0.  Doehler.  6.  Dohler.  IEEE  Transactions  on  Electron  Devices.  ED  26(10). 

1602  (1979).  - 

3)  0.  Doehler,  6.  Dohler,  International  Electron  Devices  Meeting, 
Mashlngton,  O.C.  (Oec.  1978). 
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2.0  THE  BASIC  MODEL 

The  delayed  electromagnetic  wave  in  a  slow  wave  microwave  amplifier 
propagates  in  a  space  with  an  effective  dielectric  constant  e  which  is 
determined  by  the  geometric  properties  of  the  periodic  slow  wave  structure 
and  the  space  "above"  the  line  which  contains  the  electron  beam.  In 
operation,  the  electron  beam  is  modulated  in  time  and  space  due  to  the 
electron  bunches  passing  above  the  line,  so  that  interactions  with  the  RF 
wave  can  occur.  The  volume  occupied  by  the  electron  bunches  will  therefore 
present  to  the  RF  wave  a  time  and  space  dependent  dielectric  constant.  It 
is  therefore  expected  that  the  effective  "hot"  dielectric  constant  e  of  the 
slow  wave  structure  will  also  depend  on  the  time  and  space  dependent  space 
charge  In  the  beam,  and  is  given  by: 

e  *  c  (1  +  6Ap/p0)  (2) 

Ap  *  Ap(r,t)  =  time-  and  space  dependent  space  charge 

p0  s  dc  space  charge 

6  3  geometric  factor 

If  one  assumes  now  that  the  depth  of  modulation  Ap/p0  Is  constant*,  and  that 
the  beam  is  propagating  In  z-dlrectlon,  then  Ap  can  be  developed  in  a  Fourier 
series: 

e  -  f  [  1  +  Z  6m  exp  jmu»d(t  -  z/ve)]  (3) 

ud  -  drive  frequency  (of  delayed  wave) 

ve  ■  average  electron  velocity  In  z-dlrectlon 

*  In  operation,  the  depth  of  modulation  Increases  with  z  until  saturation 
occurs.  The  assumption,  therefore,  holds  over  most  of  the  line  In  overdrive 
condition,  and  was  used  here  only  to  simplify  the  discussion. 
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The  effective  dielectric  constant  Is,  therefore,  modulated  with  all  the 
harmonics  of  the  drive  frequency.  sIf  one  simplifies  further  by  considering 
only  the  main  Fourier  component: 

e  =  r  [  1  +  6m  cos  u»p(t  -  z/ve)] 

Up  *  nwd 

then  It  Is  seen  that  the  effective  dielectric  constant  Is  modulated  with 
a  pump  frequency  fp  *  mfd.  Even  though  the  coefficient  6m  In  Relation  (4) 
may  be  small  compared  to  1,  the  existence  of  waves  at  frequencies  up  in  or 
close  to  the  stopbands.  Induces  resonances  In  each  periodic  section  of  the 
periodic  wave  structure,  which  are  more  pronounced  with  wide  stopbands,  and 
thus  may  Introduce  further  modulation  of  the  beam  at  this  frequency. 

As  required,  this  model  is  based  on  a)  the  presence  of  the  electron  beam, 
perferably  at  full  power,  b)  the  presence  of  a  drive  signal,  preferably  in 
overdrive  condition  and  c)  the  existence  of  a  periodic  structure  with  a 
stopband. 


8 


3.0  THE  CFA  DELAY  LINE 

3.1  The  Physical  Model 

Figure  4  shows  schematically  the  Interaction  space  of  a  crossed  field 
amplifier  (CFA).  The  slow  wave  structure  is  a  meander  line,  the  Individual 
fingers  of  which  are  supported  by  ceramic  bars.  The  Individual  fingers  are 
connected  to  adjacent  ones  at  opposite  ends  by  an  "empty  waveguide"  section 
(not  supported  by  ceramic  bars).  The  electron  beam  above  the  line  Is  In¬ 
jected  from  a  cathode  (not  shown),  and  remains  above  the  line  by  means  of 
a  crossed  dc  electromagnetic  field  Eq  (y-dlrectlon)  and  B0  (x-dlrection) 
which  forces  the  beam  to  move  In  z-dlrectlon  with  a  velocity  ve  =  E0/B0. 

Each  Individual  finger  of  the  slow  wave  structure  Is  therefore  perpendicular 
to  the  velocity  of  the  electrons.  As  a  result,  each  Individual  finger  n 
will  be  subjected  to  a  purely  time  dependent  dielectric  constant,  shifted 
from  the  adjacent  ones  by  a  phase  constant  phase  ±u>pP/ve,  where  P  is  the 
true  pitch  (In  z-dlrectlon)  of  the  line: 

e(n)  =  e  [  1  +  6,„  cos  up(t  -  nP/ve)] 

If  one  now  “stretches"  the  meander  line,  one  obtains  a  waveguide  (microstrip 
line)  which  Is  periodically  loaded  by  a  dielectric  of  length  £l,  as  schematically 
shown  In  Figure  5.  By  doing  this,  one  effectively  neglects  the  coupling 
between  the  bars,  but  the  frequency  position  of  the  band  gap  Is  preserved. 

Thus,  the  model  of  the  CFA,  as  shown  In  Figure  5,  Is  a  periodically  loaded 
waveguide  the  dielectric  sections  of  which  have  a  time  dependent  effective 
dielectric  constant  as  given  by  relation  (5). 


(5) 
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3.2  The  Mathematical  Model 

To  analyze  the  line  as  schematically  shown  in  Figure  5,  one  must  first 
solve  the  wave  propagation  in  a  single  section,  the  dielectric  of  which  is 
periodically  modulated  with  a  pump  frequency  fp.  Once  the  solution  has 
been  found  for  a  single  section,  boundary  conditions  at  the  interfaces  of 
the  dielectric  with  the  vacuum  section  must  be  introduced.  As  a  result,  a 
transfer  matrix  connecting  the  fields  at  the  entrance  of  a  single  dielectric 
section  n  to  the  fields  at  the  entrance  of  the  next  dielectric  section  n+1 
can  be  found.  By  multiplying  this  matrix  by  itself  N  times,  where  N  is  equal 
to  the  number  of  Individual  fingers  (approx.  96),  the  entire  transfer  matrix 
of  the  line  can  be  found. 

3.2.1  Mathematical  Solution  in  a  Dielectric  Section 

The  equations  governing  the  propagation  of  an  RF  wave  in  a  waveguide  are 
given  for  a  TEM  mode  by  Maxwell's  equations: 

3E/3x  *  -y3H/3t  3H/3x  =  +3eE/3t  (6) 

which  can  be  combined  by  Introducing  D  *  eE,  since  e  is  only  time  dependent 
In  a  single  dielectric  section 

32D/3x2  ■  ey32D/3t2 

This  equation  is  solved  by  the  method  of  variable  separation  [D  =  Dt(x)D2(t)] 
by  introducing  a  propagation  constant  6 

D"(x)/D1(x)  -  -B2  -  eyD"(t)/D2(t)  (7) 
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which  leads  to  the  solution 

Dt(x)  =  a  exp  jgx  +  b  exp  -  j&x 
together  with  Hill's  equation 

D“(t)  +  (f52/ey)D2(t)  =  0 

Using  the  dc  term  of  the  dielectric  constant  as  well  as  the  maximum  of  the 
Fourier  component  (6  «  1 ) 

e  *  e  [1  -  «  cos(u>pt  -  $)] 

one  obtains  Mathleu's  equation  by  properly  time  shifting 

D”(t)  +  82v2  (1-6  cos  wpt)Da(t)  *  0 

From  Floquet's  theorem  It  follows  that  one  solution  Is  given  by 

Fv(t)  =  eJ  tyt)  P(t)  *  ZCnejn“Pt 

where  the  characteristic  exponent  v  has  to  be  determined.  In  general,  but 
not  necessarily,  Fv(-t)  Is  also  a  solution  of  Mathleu's  equation  which  is 
linearly  Independent  of  Fv(t).  Thus  In  general 

AeJvt£CneJnwpt  +  Be-MEC^-J^  (8) 

Is  a  complete  solution  for  D2(t).  Since  |6|  Is  small  as  compared  to  one,  v 
Is  generally  real  except  for  2£v  «  nu>p,  which  will  be  excluded,  since  no  wave 
growing  exponentially  In  time  Is  observed. 
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From  the  solution  obtained  It  follows  that,  for  a  single  propagation 
constant  3(s  w/v)  an  Infinite  number  of  waves  with  frequencies  =  u  +  nu»p 
will  propagate.  The  case  u>  »  n«p/2  must  be  treated  separately. 

3.2.2  Physical  Interpretation  and  Simplification 

3.2. 2.1  Injection  of  a  Single  Wave 

Let  us  assume  now  that  a  single  wave  of  frequency  f  is  Injected  Into 
the  first  section  of  the  line.  Due  to  the  pump  frequency,  an  Infinite  number 
of  waves  with  frequencies  f n  *  f  ±  nfp  and  with  the  same  propagation  constant 
@  will  propagate  In  the  dielectric  section,  and  exit  this  section  on  both 
sides.  Consequently,  the  second  dielectric  section  will  "receive"  and  reflect 
an  Infinite  number  of  waves  with  frequencies  fn.  These  waves  will  transmit 
through  the  second  dielectric  section  with  both  the  Initial  propagation  con¬ 
stant  3  and  their  eigen-value  3n  *  2-irfn/v. 

In  principle,  therefore,  one  must  track  an  Infinite  number  of  waves  In¬ 
jected  from  both  sides  of  each  dielectric  section.  However,  since  the  amplitude 
of  the  waves  of  frequency  fn  generated  In  the  first  section,  reflected  and 
transmitted  by  this  section  decreases  rapidly  with  Increasing  value  of  n.  It 
Is  therefore  practical  to  consider  only  values  n  *  0  and  n  *  ±1. 

3. 2. 2. 2  Injection  of  Two  Waves:  Active  and  Passive  Coupling 

Let  us  now  consider  the  case  of  two  waves  with  frequencies  f,  and  f2 
Injected  Into  the  first  section.  In  this  first  section  then,  waves  with  fre¬ 
quencies  fj1  *  f x  ±  fp  and  f2*  »  f2  ±  fp  will  be  generated.  If  the  two 
primary  frequencies  f2  and  f2  are  such  that  (f,  <  f2): 

r,*t,  ■  %  (») 


14 


* 

then  these  two  waves  will  be  actively  coupled  .  If,  however 

f2  -  fi  *  fp  (10) 

then,  these  waves  are  passively  coupled.  The  differentiation  between  these 
coupling  modes  Is  made  because  two  waves 

♦l(x)eJ“lt  and  4>2(x)e^W2t 

which  are  passively  coupled  will  be  coupled  directly 
ej(wt  +  wp)t  * 

ej(“2  -  «p)t  , 

while  those  actively  coupled  will  be  coupled  Indirectly  through  their  con¬ 
jugate  complex  values 

-  wp)t  b  e-jw2t  = 

-  «p)t  =  e-jft>it  s  (ej“»tj* 

3.2. 2. 3  Presence  of  Three  Waves  and  Boundary  Conditions 
The  case  analyzed  assumes  three  primary  waves  fl#  f2,  f,  to  propagate 
In  the  CFA  delay  line,  two  of  which  are  actively  coupled  and  two  of  which  are 
passively  coupled 

fi  +  f2  -  fP  f»  -  fi  •  v  (id 

The  signification  n  ■  0,  ±1  leads  then  to  the  presence  of  six  waves  within 

*  The  terminology  of  active  and  passive  coupling  has  been  taken  from  Loulsell's 
work*).  Although  the  present  theory,  described  In  detail  In  1),  2)  and  3)  Is 
more  general  since  3  Is  not  necessarily  equal  to  /c  w/c,  this  terminology  has 
been  kept  here. 

4)  W.H.  Loulsell,  Coupled  Node  and  Parametric  Electronics,  New  York*,  Wiley  (1960). 
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_  the  CFA  lines,  three  waves  propagating  In  each  direction  (forward  and  backward). 

I  Each  of  these  waves  Is  actively  or  passively  coupled  to  other  waves,  and  Is 

reflected  or  transmitted  at  each  of  the  boundaries  between  the  dielectric  and 
the  vacuum  of  each  section.  The  mathematical  computation,  therefore,  requires 
|  the  tracking  of  six  waves,  each  wave  being  characterized  by  two  coefficients 

(amplitude  and  phase).  This  requires  12  x  12  transmission  matrices.  The 

! 

'  computer  program  used  requires  as  Input  conditions  the  geometry  of  the  delay 

line  and  the  knowledge  of  three  Incoming  waves  (amplitudes  A*(f3),  A*(f2), 

I 

A*(f  )).  Furthermore,  their  relative  phases  $lt  <j>2  and  <j>3  must  be  given, 
j  It  Is  assumed  that  the  line  Is  perfectly  matched  at  the  Input  and  output  (no 

reflections  from  the  Input  and  output  connectors).  Knowing  the  pumping  phase 
i  <j>p  required  between  two  adjacent  bars 

|  *p  *  V /ve 

the  amplitudes  A"(fx),  A“(f2),  and  A"(fs)  reflected  from  the  line  at  the  input 
as  well  as  the  amplitudes  Bj(fi)  transmitted  through  N  sections  can  be  cal- 
{  culated  (see  Figure  6). 

I  3.3  Computer  Results 

'  The  computer  results  obtained  are  summarized  herein  by  discussing  the 

|  most  pertinent  results. 

3.3.1  Results  Obtained  with  CFA  Line 

|  The  pump  frequency  fp  was  continuously  varied  between  10  and  15GHz,  since 

the  bandgap  of  the  unfolded  meander  line  was  calculated  to  be  between  12.8  and 
1  13.7GHz.  The  frequency  f,  was  varied  between  the  frequencies  200MHz  and 
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Fig.  6.  Schematic  of  boundary  conditions,  (a)  General.  (b)dj  given. 
B  7*0  set. 


fp  -  200MHz.  The  Input  conditions 
a;  -  i 

were  used,  and  the  pump  phase  <j>p  continuously  varied.  The  total  output  power 
Pi  *  (A^)2  +  (bJ)2  for  each  wave  of  frequency  fi  was  monitored.  It  became 
rapidly  clear  that  >  10  could  only  be  obtained  when  both  fx  was  close  to 
fp,  and  fp,  the  pump  frequency  In  the  bandgap,  leading  to  a  large  power  gain 
at  f2.  Varying  the  pump  phase  between  bars  lead  to  an  increase  of  P2 
above  100,  and  for  $p  =  ir,  the  result  of  Figure  7  was  obtained.  The  following 
two  conclusions  can  be  drawn: 

1)  Power  gain  is  substantial  at  low  values  of  f2(fx  -  fp). 

2)  As  the  pump  frequency  fp  Is  Increased,  the  value  f2  for  which 
P2  becomes  maximum  also  Increases 

f  m  f  _  f ,  =  f  _  f 

2max  Tp  T»max  Tp  To 
where  f0  Is  a  fixed  frequency  which  can  be  determined  from 
Figure  7  to  be 

f0  ■  12.8GHz 

which  Is  just  the  lower  edge  of  the  bandgap. 

The  Input  conditions  were  changed  by  Increasing  the  value  of  A,  from  1 
to  500,  leaving  a|  and  A*  equal  to  1.  No  appreciable  power  gain  was  obtained. 
However,  Increasing  A*  to  500,  leaving  A*  ■  A 2  -  1  led  to  a  power  gain  for  f2 
of  more  than  80dB,  see  Figure  8.  Here  again,  the  frequency  f0  ■  12.8GHz  is 
clearly  detected.  Note  also  that  about  20dB  gain  Is  obtained  at  f,  (high 
frequency). 
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Power  gain  versus  frequency  for  all  three  waves  when  fp  is  close 
to  or  in  first  bandgap. 


20 


✓ 


3.4  Experimental  Verification 

The  computer  results  obtained  predicted  that  a  very-low-frequency  com¬ 
ponent  In  the  spectrum  should  be  present  when  parametric  oscillations  are 
occurring.  The  measurement  of  frequency  components  as  low  as  a  few  megahertz 
being  impossible  at  the  output  of  the  tube  (waveguide  output),  a  coupler  was 
Inserted  at  the  Input  of  an  S-band  tube,  and  the  reflected  spectrum  measured. 
Figure  9  shows  the  low-frequency  spectrum  reflected  from  the  tube  in  over¬ 
drive  condition.  A  very  pure  signal  is  observed  at  about  50MHz  when  the  tube 
is  operated  close  to  the  ir/2  mode.  The  signal  disappears  when  the  overdrive 
condition  is  removed,  and  the  frequency  of  the  spectral  line  is  independent 
of  sole  or  cathode  voltage,  thus  refuting  the  hypotheses  that  backward  wave 
oscillations  are  responsible  for  this  low-frequency  content. 

Finally,  Figure  10  shows  the  output  spectra  under  the  same  conditions* 
(the  upper  left  picture  shows  the  spectrum  at  the  output  of  the  tube,  the 
reflected  spectrum  being  shown  in  lower  left  figure).  It  is  seen  that  the 
spacing  between  spectral  lines  is  about  60MHz  in  both  the  fundamental  and 
harmonic  spectrum,  strongly  indicating  cross  modulation  between  the  drive 
signal  and  the  harmonic  signal  with  the  low-frequency  component.  These  lines 
disappear  when  normal  drive  power  (30W)  is  applied. 

3.5  Reduced  Width  of  Bandgap  (Stopband) 

The  mathematical  model  described  above  shows  that  the  parametric  gain 
decreases  when  the  width  of  the  stopband  decreases.  It  can  be  shown  that 
excess  capacitance  or  Inductance  at  the  end  of  each  bar  of  the  meander  line 

*  The  results  shown  In  Figure  10  were  obtained  from  another  production  tube 
if  compared  to  Figure  3. 
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SPECTRUM  OF  REFLECTED  POWER 
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Fig.  9.  Spectrum  of  reflected  power. 


22 


OUTPUT  SPECTRUM 
CENTER  FREQUENCY  7J  GM* 


— 

— 

— 

— 

r — 

1 

1 

= 

- 

Ja 

d 

= 

mm 

0 

Ai 

_ 

J 

CENTER  FREQUENCY  3  9GH* 


OUTPUT  POWER  ANO  REFLECTED  SPECTRUM 
DRIVE  FREOUENCY:  39 GH* 

DRIVE  POWER:  45  Warn 

50  MHx/Ol  VISION 


REFLECTED  SPECTRUM 

Fig.  10.  Output  power  and  reflected  spectrum. 
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(finger  connections)  Increases  the  width  of  the  stopband.  Using  the  sim¬ 
plified  model  above,  the  value  6  for  the  phase  shift  per  section  can  be 
found  from: 

cos  8  =  cos  gfcj  +  ^2.  sin  b*,2 

2 

6:  phase  shift  per  section 
Jtj :  length  of  dielectric  section 

i2:  length  of  connecting  section 

y:  admittance  of  unloaded  section  &2 

Z0:  Impedance  of  loaded  section  ^ 


The  stopband  occurs  whenever  the  magnitude  of  cos  6  Is  greater  than  1.  The 
first  stopband  occurs  where  $  approaches  ir,  and  cos  8  Is  less  than  -1.  For 
larger  yZQt  the  stopband  Is  wider,  and  the  parametric  gain  will  be  greater. 

Some  experimental  cold-test  models  were  build  to  determine  the  effect 
of  possible  excess  capacitance  at  the  edges  of  the  meander  line.  The  first 
models  were  meander  lines  on  flat  continuous  substrates,  with  the  corners 
modified  as  shown  In  Figure  11.  The  effect  upon  the  width  of  the  stopband 
of  "clipping  the  comers"  to  reduce  excess  capacitance  Is  shown  In  Figure  12. 
The  point  was  never  reached  at  which  removing  more  material  would  Increase 
the  width  of  the  stopband  because  of  excess  Inductance. 


3.6  Conclusion  and  Discussion 


The  model  for  parametric  oscillations  presented  here  fulfills  the  basic 


requirements  empirically  defined  for  their  existence,  namely,  the  model  requires 
the  presence  of  an  electron  beam,  preferably  at  full  power,  and  the  presence 
of  a  drive  signal,  preferably  in  overdrive  conditions.  The  model  then  clearly 


« 


Removal  of.  material  of  meander  line  to  reduce  excess  capacitance. 


defines  a  fixed  frequency  f0  located  at  the  lower  edge  of  the  bandgap  of 
the  periodic  slow  wave  structure,  and  predicts  the  simultaneous  existence  of 
a  very-low-frequency  oscillation.  The  existence  of  this  low  frequency  os¬ 
cillation  has  been  verified  experimentally. 

The  presented  model  further  predicts  that  power  gain  Is  obtained  at 
low  frequencies  (f})  and  frequencies  close  to  the  bandgap  (f2).  This  power 
Is  obtained  from  the  pump  signal,  and  thus  apparently  from  the  second  harmonic 
of  the  drive  signal.  If  It  can  be  shown  that  power  from  the  fundamental  drive 
signal  can  be  converted  to  power  In  the  second  harmonic,  from  which  power  is 
removed,  then  the  power  "suckout"  observed  frequently  would  also  be  explained. 
Therefore,  the  necessary  further  mathematical  and  experimental  evaluation  of 
the  present  model  can  only  be  completed  If  the  direct  Influence  of  the  electron 
beam  Is  taken  Into  account. 

As  described  In  the  Interim  report,  difficulties  In  computer  programming 
occurred  during  the  first  phase  of  this  study.  It  Is  only  after  the  introduction 
of  3  coupled  waves  that  gain  could  be  obtained. 
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4.0  THE  MODEL  OF  THE  O-TYPE  TRAVELING  WAVE  TUBE 


4.1  The  Physical  Model 

Figure  13  shows  a  typical  rod  supported  helix  which  serves  as  slow  wave 
structure  In  many  0-type  traveling  wave  tubes.  The  helix,  or  spiral,  is 
supported  Inside  a  metallic  cylinder  by  a  certain  number  of  rods,  and  an 
electron  beam  of  cylindrical  symmetry  moves  with  a  velocity  ve  in  z-direction 
Inside  the  helix.  The  electromagnetic  wave.  Injected  at  one  end  of  the  helix, 
mainly  propagates  between  the  helix  and  the  cylinder  (at  ground),  and  there¬ 
fore  propagates  In  a  "quasi"  microstrip  line  with  the  velocity  nearly  equal 
to  the  velocity  of  light,  c.  To  advance  In  z-directlon  by  a  distance  P  equal 
to  the  pitch  of  the  helix,  the  RF  wave  must  travel  a  distance  l  -  /(2irR)2+P2 
along  the  helix.  Consequently,  the  wave  advances  In  z-dlrection  with  a  re¬ 
duced  velocity  v  *  cP/A,  and  synchronism  In  z-dlrectlon  between  the  beam  and 
the  wave  can  be  achieved  by  properly  dimensioning  the  helix. 

Considering  Figure  13,  It  Is  seen  that  the  wave  propagates  around  the 
helix  with  a  velocity  c  for  a  certain  distance  given  by  the  separation  between 
the  rods  supporting  the  helix.  As  a  result,  therefore,  the  RF  wave  propagates 
again  In  a  microstrip  which  Is  periodically  loaded  by  a  dielectric.  In 
Figure  13,  the  RF  wave  penetrates  and  exit1;  a  dielectric  section  three  times 
for  each  period  of  the  helix. 

Similar  to  the  case  of  the  CFA,  the  electron  beam  Is  RF  modulated  at  a 
frequency  f<|  and  contains  all  the  harmonic  frequencies  mf<j.  Following  the 
discussion  in  Section  2.0  of  this  report,  the  wave  at  point  z  will  propagate 
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Figure  13.  Schematic  of  TWT  interaction  space. 


In  a  Medium  with  an  effective  dielectric  constant  e  given  by?: 


e  ■  e  [1  +  26cos  <Dp(t  -  z/ve)]  (15) 

«p  -  mujj 

Since  the  wave  Is  slowed  down  by  the  helix  structure.  Relation  (15)  can  be 
simulated  by  assuming  a  time  and  space  dependent  dielectric  constant: 

e  =  e  [1  +  26cos  wp(t  -  z'/ve)]  (16) 

Ve  «  ve  t/P 

where  z1  Is  the  travel  distance  along  the  helix. 

Consequently,  If  one  stretches  the  helix  (as  was  done  previously  for 
the  meander  line),  the  RF  wave  propagates  In  a  microstrip  which: 

a)  Is  periodically  loaded  by  a  (small)  dielectric  section 

b)  and  whose  (effective)  dielectric  constant  varies  as  given  by 
Equation  (16). 

This  Is  schematically  shown  In  Figure  14.  Comparing  Figure  14  with  Figure  5, 

It  Is  seen  that  the  effective  dielectric  constant  of  the  helix  Is  not  only 
time  dependent,  as  It  Is  In  the  CFA,  but  Is  also  dependent  on  the  space  co¬ 
ordinate  x  (which  replaces  z').  Although  similar  In  principle,  both  models 
are  different  In  this  respect,  and  require  a  different  mathematical  model. 

4.2  The  Mathematical  Model 
4.2.1  Maxwell's  Equations 

Assuming  the  wave  to  propagate  In  the  waveguide  shown  In  Figure  14  to 
propagate  In  a  TEM  mode,  the  relevant  Maxwell's  equations  become: 

3E/3X  ■  -w  3H/3t;  3H/3x  ■  +3  cE/3t 

*  The  factor  2  in  front  of  6  was  Introduced  to  simplify  the  following  equations. 
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Microstrip  (helix) 


Dielectric  from  support  rods 


/  //// / //  ///*/'////// //l 


e  *  e[l  +  cos(ujpt-x/v  )] 


e  =  c'[l  +6'  cos  wp(t-x/Ve)] 


Figure  14.  Physical  model  of  the  helix  structure  in  0-Type  traveling 

wave  tubes  (z'  was  replaced  by  the  coordinate  x  for  simplification) 
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where  e  is  the  time  (t)  and  space  (x)  dependent  effective  dielectric  constant 
for  any  section  of  the  simulated  helix  waveguide  structure. 

For  a  given  frequency  u>,  the  solution  for  constant  e  is  obviously  given 
by  the  “free"  modes: 

exp  j  (wt  ±  8x)  8  s  u/iy  (17) 

representing  waves  with  the  propagation  constant  8  and  the  phase  velocity 
±c  ■  ±  (euj^. 

If  the  dielectric  constant  depends  on  time  and  space: 

e  *  e  [  1  +  26cos(u>pt  -  8px)]  (18) 

8p  *  V*e  =  “pp/Ave  =  {nwdP)/Uve) 

then  a  single  mode  (w,8)  will  also  contain  components  u>  +  nuip.  Using  therefore 
the  Ansatz: 


E  *  I  En  exp  j(u>  +  na»p)t  +  c.c. 

H  ■  E  Hj,  exp  j(w  +  na>p)t  +  c.c. 

where  c.c.  stands  for  the  complex  conjugate  values,  it  follows  Immediately 
by  insertion: 


dEn/dx  +  jun  yHn 
dHn/dx  +  ju,,  £En 


0 

-jwn  5e 


Ufi  *  w  +  nu>p 

j^n-1  e-J0Px  +  En+1  e^j 


(19) 


(20a) 

(20b) 


/ 
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where  the  conjugate  values  have  disappeared  since  each  of  the  equations  can 
be  modified  by  taking  Its  conjugate  complex  value.  As  discussed  previously 
for  the  CFA  delay  line.  If  a  wave  with  a  single  frequency  w  Is  Injected  Into 
the  simulated  helix  waveguide  structure,  then  waves  of  all  frequencies 
b>n  s  o>  +  nup  will  be  generated  and  will  propagate  within  the  first  section 
of  the  periodic  waveguide.  Contrary,  however,  to  the  CFA  meander  line,  these 
waves  will  not  propagate  with  a  common  propagation  constant  $.  Consequently, 
an  accurate  mathematical  model  must  keep  track  of  an  Infinite  number  of  waves 
In  both  directions.  It  is  seen,  however,  from  Equation  (20b),  that  the  waves 
generated  by  the  wave  u>  decrease  rapidly  In  amplitude  as  n  Is  increased.  It 
Is  easily  shown  that  the  amplitude  of  the  n-th  mode  related  to  oi  varies 
proportionally  to  the  n-th  power  of  6,  where  6  Is  small  as  compared  to  1 .  It 
Is  therefore  meaningful  to  consider  only  the  waves  with  n  =  0  and  n  =  ±1. 

4.2.2  Active  Coupling 

Let  us  consider  now  the  case  when  two  waves  of  frequencies  Fj  and  F2  are 
Injected  Into  the  simulated  helix  waveguide  structure.  Let  us  further  assume 
that  Fx  and  F2  are  chosen  In  such  a  way  that: 

Fi  +  F2  *  Fp 

The  Injection  of  the  wave  of  frequency  uij  will  lead  to  the  generation  of  waves 
with  frequency  Wj  -  Wp  and  u»i  +  wp.  The  value  of  u>&  -  o»p  being  equal  to  -io2, 
and  neglecting  the  term  with  frequency  wl  +  wp.  It  follows  that  the  wave  of 
frequency  Wj  Is  coupled  to  the  wave  of  frequency  w2  through  the  complex  value 
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of  the  wave  at  frequency  w2.  Similarly,  the  wave  of  frequency  ui2  will  generate 
waves  of  frequencies  w2  +  u>p  and  <d2  -  Up  s  -o>|.  Thus,  the  waves  of  frequency 
(Oj  and  o>2  are  coupled  through  their  complex  conjugate  values.  Using 


Equations  (20),  one  obtains  the  4  equations: 

dEj/dx  +  jpwiHj  =  0  (22a) 
dHj/dx  +  =  -j6eu)j(E2*  e^V)  (22b) 
dE2/dx  +  jyw2H2  »  0  (22c) 
dH2/dx  +  ja)2£E2  =  -j&D2e(Ei*  e"jBPx)  (22d) 


The  waves  (Ex,  H3)  and  (E2,  H2)  of  frequencies  wi  and  u2  are  coupled 
through  their  complex  conjugate  values.  It  can  be  shown4 )  that  the  waves  are 
actively  coupled  In  this  case,  meaning  that  power  from  (to)  the  pump  can  be 
transferred  to  (from)  the  two  waves. 

4.2.3  Passive  Coupling 

Let  us  again  assume  two  waves  of  frequency  F,  and  F3  to  be  Injected  into 
the  simulated  waveguide  structure.  Let  us  assume  now  that  the  two  frequencies 
Fj  and  F,  fulfill  the  condition: 

Fa  -  Fi  *  Fp  (23) 

The  Injection  of  the  wave  of  frequency  w3  will  lead  to  the  generation  of  waves 
with  frequency  o>3  -  top  ■  Wj  and  u>3  +  u>p.  Similarly,  the  wave  of  frequency  u>i 
will  generate  waves  of  frequency  <»>j  +  <*>p  ■  u,  and  <ui  -  u»p.  If  one  neglects 
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the  components  related  to  the  frequencies  w3  +  wp  and  u>l  -  u>p,  the  following 
relations  are  obtained  from  Equations  (20): 


dE^dx 

+ 

jywjHj 

=  0 

(24a) 

dH1/dx 

+ 

jeWjEj 

-  -jSw^E^P*) 

(24b) 

dE3/dx 

+ 

jywjHj 

*  0 

(24c) 

dH3/dx 

+ 

jui3eE3 

=  -j6eaj3(E1e“,i  BPX) 

(24d) 

It  is  seen  in  this  case  that  the  waves  (£lt  H3)  and  (E3,  H3)  of  frequencies  u»! 
and  w3  are  coupled  to  each  other  through  their  corresponding  amplitudes.  It 
can  be  shown4)  that  power  may  be  transferred  from  one  wave  to  the  other,  but 
no  power  is  delivered  by  the  pump.  Following  Loulsell4),  these  waves  are 
referred  to  as  being  passively  coupled. 

4.2.4  Relationship  to  Normal  Modes 

Relations  (24a)  and  (24b)  can  be  combined  by  multiplying  the  second  relation 
by  the  impedance  Z0  ■  Ai/l  of  the  waveguide.  If  one  adds  and  subtracts  the 
resulting  relation  from  Equation  (24a),  one  obtains  the  relation: 

*i±  ♦»*  =  -j  f  (4>3+  +  ♦,  J  ejBPx  (25) 


where  we  have  introduced  the  normal  modes: 


and  the  differential  operator 


J 


V  ■  fc  *  * 

Equation  (26)  is  the  classical  representation  of  normal  modes.  The  power 
carried  by  the  mode  Is  |<J>  |2  -  |i|>n_|2,  Indicating  that  $n+  Is  a  forward 

mode,  and  <|>n_  is  a  backward  mode.  This  can  be  easily  seen  from  relation  (25) 
when  neglecting  coupling.  One  obtains  then: 

♦  l±  -  <=±JS'X 

which  corresponds  to  waves  with  positive  and  negative  phase  velocity. 

4.2.5  Parametric  Amplification 

Let  us  assume  two  waves  of  frequencies  Wj  and  o>2  to  propagate  In  the 
waveguide,  and  let  us  assume  that  they  are  actively  coupled: 

“l  +  U2  =  “p 

and  that  their  uncoupled  propagation  constants  are  Bx  and  82»  such  that: 

Bp  =  Bt  +  B2  +  8 

The  corresponding  normalized  coupled  mode  equations  are  then  given  by: 

A1+  <J>i+  -  -j|  /ev  <fr2*  e-j6Px 

A2-  <J>2+  -  J  g  ^  “a  e^pX 


(27) 


(28) 
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where  the  conjugate  complex  of  the  second  equation  was  taken  (A+*  »  A_). 
In  normal  mode  nomenclature,  the  coupling  coefficient  c  from  wave  1  to 
wave  2  and  the  coupling  c2i  from  wave  2  to  wave  1  are  given  by: 

c12  =  -j  |  uj 

c2i  =  j  J  “2 

showing  immediately  that  the  Manley  Rowe  equations  are  fulfilled: 

Cn  C21* 

0)j  0)2 

Introducing: 

Bi'  *  Ba  +  38 

B2'  *  P2  +  jb 

and  using  the  Ansatz: 

4>j+  *  Aj+  e“JBj'x  j  ■  1.  2 

one  obtains  by  Inserting  Into  the  system  (28): 


(ax  "  4  ) 


’i+ 


Ci2A, 


2+ 


(fc*4) 


'2+ 


C21*l 
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Solutions  of  this  system  of  equations  Is  easily  obtained  by  assuming  A1+ 
and  A2+  to  be  proportional  to  exp  Xx.  The  Eigen-value  X  must  satisfy  the 
relation: 

X  *  ±  j  ^  ~6%Z~2  -  0*  (29) 

showing  that  exponential  growing  waves  are  obtained  for  very  small  values  of  0. 

The  energy  gained  by  the  waves,  therefore,  must  be  obtained  from  the  pump. 

For  0  =  0,  the  maximum  gain  Is  obtained,  namely  Xm  =  I  6/0 102. 

It  Is  Interesting  to  note  that  exponential  gain  with  a  backward  wave 

can  be  obtained.  For  0  *  -20z  namely,  relation  (29)  can  give  exponential 

gain: 

X  -  ±  B2  /6*0 j/02  -  4 

However,  the  maximum  gain  becomes  very  small. 

4.2.6  Frequency  Conversion 

Let  us  now  consider  the  case  of  two  waves  of  frequency  F,  and  F,  which  are 
passively  coupled: 


F,  -  F,  -  Fp  (30) 

Let  us  further  assume  that  these  waves  have  uncoupled  propagation  constant  0, 
and  03  such  that: 

6,  -  Bi  ♦  8  »  Bp 


The  corresponding  normalized  coupled  mode  equations  are  easily  obtained  from 
equation  (24),  namely: 


^1+^1+  s  ” j  2  W1  $3+ 

^3+^3+  a  -j  f  J&  w3  4>1+  e"^BPx 

As  previously  shown,  the  Manley  Rowe  equations  are  fulfilled.  This  system  of 
equations  can  be  solved  by  using  the  approach  described  In  Section  4.2.5. 

We  will  employ  Laplace's  transform  to  demonstrate  a  more  general  technique  to 
solve  such  a  system  of  equations. 

If  ^j+(p)  Is  the  image  of  4>j+(x),  It  follows: 

(p  +  jBi)  $i+(p)  *  »iC  $,+(p  -  jBp) 

(p  +  JB3)  4s+(p)  *  «3c  <h+(P  +  JBp) 

From  the  second  of  these  relations,  the  value  $3+(p  -  jBp)  can  be  found 
by  replacing  p  by  p  -  j&p.  It  follows  then  that  $3  can  be  eliminated,  leading 
to  an  equation  containing  $1+  alone.  This  relation  can  only  be  fulfilled  If: 

P2  +  Jp(&i  +  B,  -  Bp)  -  Bj(B3  -  Bp)  -  c2u)x«,  *  0 

leading  to: 

P  B  -JtBi  +  B,  -  Bp  ±  /(Bj  -  B,  +  Bp)2  -  4c2WjWS]  /  2 

Since  c2  Is  negative*.  It  follows  that  p  Is  always  Imaginary,  so  that  no  In¬ 

stabilities  (exponential  growth)  can  occur.  In  this  case,  power  Is  transferred 
back  and  forth  from  the  two  distinct  modes  $l+  and  4S+. 

*  c  Is  purely  Imaginary 
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4.2.7  Boundary  Conditions 

Considering  the  periodically  loaded  waveguide  shown  In  Figure  14,  It  Is 
seen  that  the  waves  alternately  will  propagate  In  an  empty  section  of  the  wave¬ 
guide  followed  by  a  ''small"  section  filled  with  a  dielectric.  At  each  of 
the  Interfaces  between  the  empty  waveguide  and  the  filled  waveguide,  boundary 
conditions  must  be  fulfilled.  These  boundary  conditions  require  the  matching 
of  the  field  components  {E,  H)  for  any  given  frequency.  From  the  definition  (26) 
of  the  normal  modes.  It  follows  that  the  two  functions: 

Ej  *  2  ^ o  (4j+  +  4>j.) 

Hj  3  2(4»j+  - 

must  be  continuous  at  each  of  the  boundaries  considered.  Note  that  j  Indexes 
the  frequency  considered. 

4.2.8  Conclusion 

We  have  established  the  mathematical  fundamentals  of  the  model  proposed  for 
the  traveling  wave  tube  In  terms  of  the  coupled  mode  theory.  As  In  the  case 
of  the  CFA,  we  will  consider  waves  with  three  different  frequencies,  which 
means  that  we  will  have  to  track  6  waves,  3  In  each  direction.  The  boundary 
conditions  (31)  must  then  apply  at  each  of  the  many  Interfaces  between  the 
empty  and  filled  waveguide  sections. 


(31) 


4.3  Three  Waves  ‘ 

4.3.1  Basic  Relations 

We  shall  consider  three  waves  of  frequency  Fx,  F2  and  Fs  to  propagate 
In  each  direction.  We  shall  assume  that  the  waves  1  and  2  are  actively 
coupled,  and  that  the  waves  1  and  3  are  passively  coupled,  l.e. : 

F,  +  F2  =  Fp  (32) 

F3  -  F,  -  Fp 

In  view  of  the  results  obtained  In  Section  4.2  of  this  report,  it  Is  seen 
that  the  "free"  propagation  constants  6l  *  w-f’^u  fulfill  equations  similar 
to  Equation  (32),  namely: 

3j  +  &2  a  *  6p  (33) 

B3  -  B,  *  o)p*/eu  =  £p 

where  we  have  Introduced  which  Is  not  necessarily  equal  to  Bp,  the  prop¬ 
agation  constant  of  the  "dielectric  wave",  which  propagates  with  a  velocity 
v  determined  mainly  by  the  velocity  of  the  electrons. 

Considering  the  results  obtained  In  Section  4.2,  It  Is  seen  that  a  system 
of  6  differential  equations  of  the  first  degree  must  be  solved: 

A,±  4>1±  -  *jcBi(<fr2+*  +  4>2.*)  exp  -jBpX  (34a) 

+jc8,(«,+  +  $,J  exp  JBpX 
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A2?  *2±*  =  ±  JC®2  ($1  +  +  <(>1-)  eXP  j&pX 

a9±  <frs±  s  +  jc6s  ($1+  +  4>,J  exp  -jBpX 

A-j±  =  (d/dx)  ±  jB-ji  c  *  6/2 


The  first  step  In  analyzing  the  proposed  model  will  be  to  solve  this  system 
of  equations. 

4.3.2  The  Dispersion  Equation 

The  system  of  equations  (34)  can  be  Laplace  transformed: 


♦(p)  *  /  e“px  $(x)  dx 

o 


Six  equations  are  then  obtained,  which  contain  the  Image  functions  4>j±,  some 
of  which  with  arguments  p  ±  jBp.  These  Image  functions  can  easily  be  elim¬ 
inated  by  shifting  p  to  p  ±  jBp.  Eliminating  the  functions  $2±  and  $3± 
leads  to  two  equations  for  $1±: 


-IC 


*1+  (P) 


P  -  j0,  +  K 


♦l-  (P) 


0 


which  can  only  be  fulfilled  If  the  determinant  of  the  2x2  matrix  vanishes. 
This  leads  to  particular  values  of  p  ■  -jB,  which  must  fulfill  the  dispersion 
equation: 


b2  -  b,2  »  «aBi2 


02* 


(B-Bp)2  -  Ba*  (B+Bp)2  -  B,s 


42 

/ 


(34b) 


(35) 


If  6  =  0,  the  "free"  solutions  6  ±  Bx  are  obtained.  If  5  f  0,  then  relation  (35) 
represents  a  6th  order  equation  leading  to  6  different  solutions.  This  could 
be  predicted  since  we  started  with  a  system  of  6  differential  equations  (34). 

In  general*  62  is  very  small  as  compared  to  1.  Unless  at  least  one  term 
of  the  denominators  of  the  right  hand  In  (35)  nearly  vanishes,  the  solutions 
of  equation  (35)  are  given  approximately  by  B  *  ±BX.  "Strong"  deviations  from 
these  solutions  can  only  occur  If  and  only  If  one  or  both  of  the  following 
relations  is  satisfied: 

±6,  ±  B2  »  Bp 

±8X  ±  Bj  =  Bp 

Remembering  that  Bj  are  positive  values,  and  considering  relation  (33),  it 
follows  that  Interaction  occurs  when  either  or  both  the  following  relations 
are  satisfied: 


Pi  +  &2  s  6p  *  &p 

Bj  "  B1  *  Bp  *  Bp 

The  right  hand  side  of  equation  (35)  will  be  small  as  compared  to  1  (pro¬ 
portional  to  62)  If  B  Is  set  approximately  equal  to  -Bi-  Consequently,  one 
can  deduce  that  the  backward  wave  1  will  never  be  "strongly”  coupled  to  the 
other  waves,  and,  as  such,  can  be  assumed  to  propagate  "free"  with  the  un¬ 
perturbed  propagation  constant  -Bx.  Similarly,  one  can  deduce  that  the 


I 
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backward  waves  2  and  3  will  also  propagate  "freely**  with  their  respective 
propagation  constants  -02  and  -0,  respectively.  Thus,  one  can  assume  the 
backward  waves  to  propagate  as: 

<J»1_  ~  exp  +j0ix 

If  one  neglects  now  the  respective  coupling  terms  related  to  backward  waves 
In  the  system  of  equations  (34),  one  finds  a  dispersion  equation  of  the  third 
order  only,  namely: 

03  _  S2 

A0  +  A0p  A0  -  A0p 

A0  -  0  -  0ji  A0p  =  0p  -  0p 

This  relation  merits  some  discussion.  Consider  the  following  cases: 

(a)  Neglect  wave  3  (02  ■  0) 

The  solution  corresponds  then  to  the  solution  obtained  earlier  for 
actively  coupled  waves  1  and  2: 

A0  «  (A0p/2)  ±  j  V620,02  .  ABp2  1  2 

Under  synchronism  (A0p  ■  0),  maximum  gain  for  the  exponential  waves  are  obtained. 

(b)  Neglect  wave  2  (02  ■  0) 

The  solution  corresponds  to  the  solution  obtained  earlier  for  pas¬ 
sively  coupled  waves  1  and  3: 

A0  -  -(A0p/2)  ±  j  40p*  +  «20,0,  /  2 

In  this  case,  no  growing  waves  are  obtained. 


(36) 
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(c)  Synchronism  (gp  ■  Bp) 

In  this  case.  It  Is  easily  seen  that  the  solutions  become: 

Bj  and  Bj  ±  s;  s  »  (6/2)*'B1(Bs  -  B2) 

Since  It  can  be  shown  that  Bs  -  B2  *  2BX,  no  growing  waves  can  be  expected. 

In  view  of  the  growing  waves  one  obtains  between  actively  coupled  waves, 
this  result  seems  surprising.  However,  the  passively  coupled  mode  dominates 
the  actively  coupled  mode,  and  power  Is  essentially  converted  among  the  various 
waves. 

(d)  General  case  (assynchronlsm) 

In  the  general  case  of  assnychronlsm  (ABp  t  0),  It  can  be  shown  that 
growing  waves  will  exist  If  the  synchronism  factor  ABp  Is  slightly  positive 
or  negative.  Under  certain  conditions,  discussed  In  Appendix  1,  the  Imaginary 
part  of  B  becomes: 

Im(B)  *  ±  .24  6Bq 

(6q/B,  »  1) 

In  general,  however,  one  must  solve  the  dispersion  equation  (36)  to  obtain 
the  exact  values  of  the  solutions  B^< 


4.3.3  The  General* Solution 


Representing  the  solutions  of  the  dispersion  equation  (36)  to  be 


A0  —  S  j »  S2  >  Sj 


then  the  most  general  solution  for  the  three  waves  Is  obtained  by  using  the 
ansatz: 


$1+ 


3 

I 


An  e"^SnX 


e'JM 


(38a) 


Inserting  this  value  Into  relations  (34b)  and  (34c)  and  keeping  In  mind  that 
coupling  from  and  to  backward  waves  Is  neglected,  leads  to: 

-jsnx 


$ 


2+ 


*  ee2 


Z  An  e 
1 


exp  j(02  +  ABp)x 


(38b) 


A8p  -  *n 


$3+ 


C8S 


Ar  e~jSnX 
A6p  +  sn 


exp  -j(B,  +  ABp)x 


(38c) 


where  the  coefficients  An  are  determined  by  boundary  conditions.  For  example, 
the  knowledge  of  4>i+(0)  for  all  three  waves  (1  ■  1 ,  2,  3)  at  x  *  0  uniquely 
determines  An. 

4.3.4  Conclusion 

The  most  general  solution  has  been  determined  when  3  waves  of  frequencies 
Fj,  F2  and  Fs  are  propagating  In  each  direction  In  a  given  section  of  a  wave¬ 
guide.  Two  of  these  waves  were  assumed  to  be  actively  coupled,  and  two  of 
these  waves  were  assumed  to  be  passively  coupled.  The  simultaneous  presence 
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of  all  three  waves  does  not  lead  to  exponential  growing  waves  when  perfect 
synchronism  is  assumed.  A  certain  value  of  the  asynchronism  factor  A6p,  how¬ 
ever,  leads  to  exponential  growing  waves.  The  sensitivity  of  the  synchronism 
factor  ABp  on  the  exponential  growth  of  the  waves  may  explain  the  extreme 
sensitivity  of  the  "power  suck-out"  to  external  and  internal  tube  parameters, 
such  as  beam  voltage  and  assymmetry  in  the  placement  of  the  helix  supporting 
rods. 

A  generalized  theory  with  these  three  waves  requires  solving  a  dispersion 
equation  of  the  six-th  order.  Simplification  results  when  no  coupling  from 
and  to  the  backward  waves  is  assumed.  The  third  order  dispersion  equation  leads 
to  three  solutions,  two  of  which  can  be  complex  conjugate.  The  general  solution 
(38)  must  then  be  applied  to  each  section  of  the  waveguide,  and  proper  boundary 
condition  introduced.  This  requires  again  a  complex  computer  code,  and  the 
results  will  not  be  explicative  from  the  physics  point  of  view.  Since  two 
actively  coupled  waves  lead  to  exponential  growing  waves,  we  will  apply  In 
the  next  section  the  boundary  conditions  to  waves  Fx  and  F2  which  are  actively 
coupled  (Section  5.0). 


5.0  PARAMETRIC  TWT  AMPLIFICATION 


5.1  General 

In  order  to  understand  the  physics  involved  in  both  the  cases  of  CFA  and 
the  TWT,  we  will  consider  in  this  section  two  actively  coupled  waves  such  that: 


Uj  +  (»2  -  (Dp 

6i  +  e2  55  0P 

It  follows  then  from  the  discussion  In  Section  4.0  that  the  general  solution 
of  the  system  containing  4  waves  (2  in  each  direction)  Is  given  by: 

4>i+  =  (aj  esx  +  a2  e”sx)  exp  -j&jX 


bx  exp  j8iX 

jc^|  (a1  esx  -  a2  e’sx)  exp  +iB2x 
b2  exp  -jp2x 
s  ■  c  /SjBj  *  | 


(40a) 

(40b) 

(40c) 

(40d) 


where  a^  and  b^  are  constants  of  Integration. 

' ?t  us  define  the  vector: 

*  -  (<fr1+.  <>!_.  42*»  ♦.!>  (41) 

It  is  then  easy  to  define  the  transformation  matrices  for  a  single  section 
of  the  waveguide. 
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5.2  The  Transformation  Matrices 


5.2.1  Transformation  Through  a  Single  Empty  Waveguide  Section  of  Length  fti 
Figure  14  shows  the  schematic  of  the  waveguide  to  be  analyzed.  Let  us 

assume  that  the  "entrance”  of  a  single  empty  waveguide  section  Is  located  at 
x  =  0,  so  that  the  "exit"  of  the  empty  waveguide  section  Is  located  at  x  « 
see  Figure  15. 

If  the  value  of  the  vector  ij)  at  x  =  0  is  known  (see  definition  41),  then 
the  value  of  i|;(£i)  in  the  empty  waveguide  can  be  easily  determined  by  using  the 
relations  (40).  One  obtains  a  transformation  matrix  (Li)  which  transforms  the 
vector  such  that: 

JUi)  8  UiMO)  (42) 

where  Lj  is  defined  in  Figure  16.  This  matrix  can  be  decomposed  into  a  matrix 
(L j 0)  representing  the  transmission  In  the  absence  of  pumping  (s  =  0)  and  a 
coupling  matrix  (5X)  which  describes  the  coupling  through  the  pump  frequency  up. 

5.2.2  Interface  Boundary  Condition 

The  boundary  conditions  (31)  must  now  be  applied  at  the  Interface  between 
the  vacuum  and  the  dielectric  (see  Figure  15).  The  vectors  at  the  position 
x  =  -  A  and  at  the  position  x  ■  l-i  -  A  (A  -*•  e)  can  be  related  to  each  other 

by  a  transformation  matrlc  (B0): 

$(*i  -  A)  =  (B0)  *  (ij  +  A)  (43) 

where  (B0)  is  defined  in  Figure  17.  We  are  now  In  the  position  to  relate  directly 
^(tj)  In  the  dielectric  to  the  Initial  vector  ifi(0)  In  the  vacuum  by  the  trans¬ 
formation: 

J(*x)  »  (B0-‘)  (Lt)  ?  (0) 

Note  that  B0  Is  a  transformation  matrix  which  does  not  depend  on  the  pump  wave  Up. 
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Figure  16.  The  transformation  matrix  (L,)  which  can  be  represented  by  the 
product  (L,)  •  (Lj,)(s,). 
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Figure  17.  The  4x4  boundary  transformation  matrix  (B0)  decomposed 
Into  2x2  matrices. 
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5.2.3  Transformation  Through  the  Dielectric  Section 
The  behavior  of  the  solution  In  the  dielectric  portion  Is  similar  to  the 
behavior  of  the  solutions  In  the  vacuum  section  If  one  replaces  6^  by  ki, 
where  *  Bi/ZeJ.  being  the  dielectric  constant  of  the  dielectric  material. 

To  be  exact,  the  solution  of  the  dispersion  equations  will  be  ±S,,  where  s, 

Is  either  real  or  Imaginary,  depending  on  the  synchronization  factor.  There¬ 
fore,  the  transformation  matrix  through  the  dielectric  section  (L2)  will  be 
similar  to  the  transformation  matrix  (Lt ) .  We  will  assume  that  12  Is  so 
small,  that  the  coupling  matrix  (S2)  Is  nearly  equal  to  unity: 

$(*,  +  *2)  a  U2)  *  (tx)  (44) 

(L2)  *  (Lao)  (S2)  «  (L20) 

where  the  transformation  matrix  (L20)  Is  shown  In  Figure  18  (and  should  be 
compared  to  Figure  16). 

5.2.4  The  Transformation  Through  a  Complete  Single  Section 
If  one  applies  a  boundary  condition  transformation  similar  to  the  one 
given  by  relation  (43),  one  finds  that  the  transformation  one  has  to  apply 
Is  just  given  by  Bo.  Consequently,  It  follows  that  4»(Jtt  +  12)  *  $(p)  at  the 
entrance  of  a  vacuum  section  can  be  related  to  the  vector  |(0)  of  the  fol¬ 
lowing  section  through  a  transformation  matrix  given  by: 

&P)  «  (Bo)  (L20)  (S2)  (B0‘»)  (Lio)  (SJ  J  (0)  (45) 
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Figure  18.  The  4  x  4  transmission  matrix  (L2)  ■  (L2B)  (s,)  «  (L,.) 
through  the  dielectric. 
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which,  as  discussed  earlier,  can  be  simplified  by  setting  (S2)  *  (I): 
ftp)  ■  (T0)  (S,)  fto) 

(T0)  •  (B0)  (Li.)  (Bo*1)  (L„) 

where  (T0)  is  the  unperturbed  transformation  matrix  one  obtains  when  coupling 
(pumping)  between  the  waves  Is  neglected. 

5.2.5  The  Complete  Transformation  Matrix 

It  Is  clear  that  If  N  Individual  sections  are  present,  then  one  can 
relate  the  Input  vector  iji(0)  to  the  helix  to  the  output  vector  $(Np)  of 
the  helix  by  successive  transformations  (46)  to  obtain: 

ftNp)  •  (T„S1)Nft0) 

We  note  here  that  (T0)  Is  the  Inperturbed  transformation  matrix,  and  that 
(T0)  Is  multiplied  by  a  coupling  matrix  (Si)  which  Is  nearly  equal  to  the 
unit  matrix  stx  «  1.  It  Is  therefore  of  Interest  to  discuss  the  matrix  To 
before  solving  relation  (47). 

5.3  The  Unperturbed  Matrices 

In  the  absence  of  the  coupling,  l.e. ,  (S,)«I,  both  waves  of  frequencies  u>j 
and  o>2  will  propagate  uncoupled  through  the  periodic  waveguide  structure.  Con¬ 
sequently,  the  4  x  4  matrix  (T0)  can  be  reduced  to  2  x  2  matrices  Z(i)  such  that: 

Z(l)  0 

To  * 

0  2(2)/ 


(47) 


(48) 
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where  Z(1)  represents  the  transformation  of  a  single  wave  of  frequency 
The  matrix  Z  can  be  calculated  rather  easily,  and  one  finds  the  matrix  shown 
In  Figure  19.  This  matrix  fulfills  the  condition  Det  Z  ■  |Z|  ■  1,  and,  fol¬ 
lowing  Cowley  Hamilton's  relation,  can  be  represented  by: 


Z  ■  I  cos  0  ♦  K  sin  0 


where  0  Is  defined  by: 


cos  0  *  ^  (Zu  +  Z22)  *  Real  Zn 


cos  @4,  cos  k4,  -  L-t-§.  sin  04,  sin  k4, 
2/e 


and  K  given  by: 


K  sin  0 


Zn  -  Zj2 
2 


-  Zii-z22 


It  can  be  shown  rather  easily  that: 


K2  ■  -  I;  Zn  ■  I  cos  n0  +  K  sin  n0 


where  I  Is  the  unit  matrix. 

Relation  (51)  describes  the  w-0  diagram.  If  the  right  hand  Side  of 
Equation  (50)  Is  larger  than  or  equal  to  1,  then  0  becomes  Imaginary,  and  no 
propagation  can  occur  In  a  waveguide  made  up  of  many  Individual  periodic  sections. 
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If  6  Is  real,  then- relation  (51)  describes  Floquet's  theorem:  the  waves 
suffer  a  net  phase  shift  6  per  single  section. 

Of  Interest  for  later  use  are  the  eigenvalues  X  and  the  eigenvectors  t 
of  the  matrix  Z.  It  Is  easily  seen  that  the  relation: 

zfn  *  *n  fn 
can  be  verified  for: 

*1  *  X2*  -  ej6  (53) 

I1>2  =  (-*12,  Kxl+j) 

showing  that  each  wave  E1>2  Is  phase  shifted  through  the  N  sections  of  the 
helix: 

t  .  .  ±jn6  1 

5i  2  (out)  «  e  Cx  *  (In) 

t  * 

Note  that  the  eigenvalues  fulfill  the  second  order  equation: 

X*  -  2X  cos  0  +  1  =  0  (54) 


It  follows  therefore  that  the  unperturbed  matrix  T0  In  equation  (48)  has 
the  4  eigenvalues  Xl#  X2  ■  Xj*.  X,,  X„  *  X,*,  where  X1>2  are  the  eigenvalues 
associated  with  the  matrix  Z(1),  l.e. ,  the  frequency  u^,  and  X,  „  are  the 
eigenvalues  associated  with  the  matrix  Z(2),  l.e.,  the  frequency  w2.  The 
eigenvectors  of  the  matrix  T0  are  then  given  by: 


Cx, a  »  [-  K12(l),  K,i(l)*j,  0,  0] 

Clf,  -  [0,  0,  -K12(2),  Kn(l)*j] 

X  -  eij01-  X  -  e*j01 

”1,8  *  »  A*,*  * 


(55) 
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|  5.4  General  Method  of  Solution 

"  The  general  method  which  will  be  applied  to  solve  relation  (47)  will  be 

^  as  follows: 

1)  determine  the  4  eigenvalues  of  the  matrix  T  *  T0S1 

2)  determine  the  4  eigenvectors  o|  of  the  matrix  T  » 

3)  decompose  the  Input  vecotr  $(0)  Into  4  eigenvectors  E  on 

4)  apply  the  multiplication  rule: 

Tn  *(0)  «  Z  %  aJ[  on 

5)  compare  the  resulting  vector  to  the  vector  $(Np) 

When  applying  the  multiplication  rule,  let  us  assume  that  Ar  *  Xn(l  +  e„), 

then: 

TN  *p(0)  «  I  c*  xj  (1  +  e*)"  on 

so  that  gain  can  be  expected  whenever  |1  +  en|  >  1,  since: 

(1  +  e)N  «  exp  [N  £n  (1  +  en)]  *  exp  Ncn 

'One  can  therefore  expect  a  power  gain  G  whenever  the  real  part  of  en  Is  positive. 

The  gain  G  will  then  be  given  by: 

Gn  «  8.7  N  Real  e,,  (56) 
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5.5  Results  and  Discussion 

It  can  be  shown  that  the  4  eigenvalues  Ar  of  the  transformation  matrix  must 
fulfil  a  4th  order  equation  which  Is  discussed  In  Appendix  II.  This  4th  order 
equation  with  complex  coefficients  has  been  numerically  solved  by  computer. 

With  N  -  100,  small  gain  of  a  few  hundreds  to  a  few  tenths  of  a  dB  were  generally 
obtained*  except  In  the  vicinity  of  61(2  ~  kir.  In  this  case,  a  double  or  quadruple 
solution  was  expected,  and  Inaccuracies  In  the  computer  code  were  evident.  How¬ 
ever,  gains  in  excess  of  100  dB  were  calculated. 

A  general  discussion  on  the  eigenvalues  Is  made  In  Appendix  II.  Of 
primary  importance  are  the  two  parameters  et  and  e2,  which  describe  the  actual 
phase  shift  of  the  waves  Fi  and  F2  through  a  single  waveguide  section.  Dif¬ 
ferent  cases  are  considered. 

Case  I  -  General  Case:  kir  f  6i  f  82  f  ntr 

The  gain  obtained  In  this  general  case  Is  relatively  small,  being  pro¬ 
portional  to  the  square  of  the  amplitude  of  the  dielectric  modulation  (l.e.  ~  62). 
Since  6  Is  small  compared  to  1,  no  appreciable  gain  Is  obtained.  This  was 
verified  by  computer  code. 

Case  II  -  The  it  modes  (9i  or  8a  °  kir,  but  9i  i  8,) 

It  Is  shown  In  Appendix  II  that  If  one  of  the  waves  Is  a  it  mode,  then  a 
gain  proportional  to  6  Is  obtained.  This  leads  to  a  gain  of  a  few  dB's,  as 
verified  by  the  computer. 

Case  III  -  Dual  and  Equal  w  modes  (kir  ■  9i  ■  9,  +  2nir) 

Maximum  gain  Is  obtained  when  both  waves  are  ir  (or  nir)  modes.  The  gain  here 
Is  found  to  be  proportional  to  the  square  root  of  S  and,  with  the  values  6  •  10"2, 
e  “  5,  W&i  “  10%,  gains  In  excess  of  100  dB's  were  calculated. 

*  Computer  evaluation  was  made  by  using  6  «  10"2,  e  «  5,  laft\  »  10% 
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As  shown  In  Appendix  II,  the  relevant  values  of  en  *  Xn(l  +  en)] 
are  approximately  given  by: 

en2  s  ±  slnhfd/SiSa  !i/i2)A 

2/e 

A2  a  (slng^!  s1ne2ti) 


Relation  (57)  has  been  confirmed  within  20%  error  by  computer  calculations. 

tf2 

It  Is  surprising  to  find  the  gain  to  be  proportional  to  61 

Case  IV  -  Dual  and  Opposite  w  inodes:  (kvr  *  0!  *  02  +  (2n  +  1  )ir) 

The  gain  one  obtains  in  this  case  Is  porportlonal  to  6,  and  corresponds 
to  the  single  ir  mode  case.  This  case  Is  Important,  however,  since  It  can  directly 
be  related  to  the  results  obtained  for  the  CFA. 

In  the  case  of  the  CFA,  maximum  gain  was  obtained  when  F,  was  first  at 
the  lower  edge  of  the  first  band  gap  (0,  *  ir).  A  very  low  frequency  oscillation 
(F2  =  50  MHz)  was  then  determined.  This  corresponds  to  02  *  0. 

It  can  therefore  be  expected  that  even  higher  gains  can  be  obtained  In 
the  case  of  the  CFA  Is  both  6j  and  02  are  set  equal  to  kvr  modes  (e2  *  kir, 
e2  *  e,  +  2nir).  This  case  was  not  analyzed  as  yet:  the  approximations  made  for 
the  CFA  evaluation  do  not  allow  to  set  Ft  ■  F2. 


6.0  CONCLUSIONS  AND  FUTURE  WORK 

The  general  model,  proposed  by  Northrop  and  supported  by  AFOSR,  for  the 
occurence  of  parametric  oscillations  In  high  power  microwave  amplifiers  has 
been  applied  to  both  Injected  beam  crossed  field  amplifiers  and  traveling 
wave  tubes.  The  results  obtained  clearly  Indicate  that  the  model  can  describe 
the  effects  of  one  class  of  spurious  oscillations  observed  experimentally. 

In  the  case  of  the  CFA,  a  constant  frequency  located  within  the  first 
band  gap  has  been  Identified.  The  occurence  of  the  spurious  oscillation  has 
been  linked  to  a  very  low  frequency  oscillation.  This  low  frequency  os¬ 
cillations  has  been  observed  experimentally,  and  was  shown  not  to  be  correlated 
to  backward  wave  oscillation.  Experimentally  the  band  gap  of  meander  lines 
has  been  reduced  by  properly  matching  the  finger  end  connections.  The  results 
have  been  published  In  2  publications2*3. 

The  model  of  the  TWT  was  simplified  In  order  to  obtain  a  better  under¬ 
standing  of  these  oscillations.  The  occurrence  of  simultaneous  it  modes 
leads  to  gains  of  the  order  of  100  dB's  or  more.  The  analysis  carried  out  In 
this  report  Impllcltely  assumes  perfect  synchronism  (Bi  +  B2  *  8p).  In 
general,  this  Is  not  the  case.  Appendix  III  shows  the  dispersion  equation  for 
the  assynchronlsm.  A  detailed  analysis  should  be  made  ustng  the  proper  dis¬ 
persion  equation, 

A  basis  for  the  occurence  of  parametric  oscillations  has  been  .established. 
If  It  can  be  shown  that  power  from  the  fundamental  drive  signal  can  be  converted 
to  power  In  the  second  or  higher  order  harmonic  from  which  power  Is  removed  to 
build  up  the  parametric  waves,  then  the  "power  suck-out"  observed  frequently 
would  be  explained.  This  subject  should  be  Investigated,  together  with  a 
correlation  between  experimental  data  and  the  theoretical  predictions  made  In 
this  report. 
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APPENDIX  1 

The  Dispersion  Equation 


The  dispersion  equation  (38)  can  be  written  In  normalized  form: 

x3  +  3px  +  2q  =  0  (IJ) 

p  =  -y  [y2  +  l];  q  -  yz 

x  =  );  y*  V2A&P/6B!-,  z  -  Bq/26, 

The  relevant  parameter  to  discuss  is: 

*  ■  p3  +  q2  ■  -27  (y2  +  i)3  +  y2z2  (1.2) 

If  k  is  positive,  then  one  real  and  2  conjugate  complex  solutions  exist. 

Let  us  therefore  discuss  the  possibility  of  conjugate  complex  solutions. 

In  order  to  obtain  growing  waves,  one  must  fulfill  the  condition: 

<  y2z2  (1.3) 

In  order  to  discuss  this  condition,  let  us  discuss  the  value  for  k.  For  y  *  0, 
k  Is  equal  to  -1/27,  and  Is  negative.  The  maxima  (and  minima)  of  k  as  a  function 
of  y  are  determined  by  the  condition: 

3k  _ 

—  ■  0 
3y 


which  leads  to  the  following  conditions: 
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(a)  y  -  0  k  <  0 

(b)  y2  +  1  =  ±3z 

Condition  (a)  Is  without  direct  Interest,  since  k  <  0.  Condition  (b)  requires 
y2  to  be  positive,  and  therefore: 

y2  =  3z  -  1  >  0 

Inserting  this  value  for  y  Into  (1.2)  leads  to  the  condition: 

=  *2  (2z  -  1) 

which  is  always  fulfilled  since  2z  *  Bq/Bi  and  the  waves  are  actively  coupled. 

It  follows  therefore  that  exponential  growing  waves  will  always  exist  If  the 
"synchronism"  factor  y  is  set  equal  to 

y  *  ±  ^3z  -  1 

and  It  Is  clear  that  a  certain  "bandwidth"  A Bp  exists  around  this  value 
in  which  growing  waves  exist.  If  k  Is  positive,  then  one  usually  defines 
the  two  parameters  u  and  v  given  by: 

u*  *  -q  +  /c;  v*  *  -q  -  (1.4) 

and  the  3  solutions  are  then  given  by: 

u  +  v  .*^3#  * 

■  u  +  v;  x2fl  - - —  ±j~2<u  -  v) 
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(1.5) 
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Maximum  gain  Is  therefore  obtained  when  k  Is  maximally  positive,  1 
k  =  <m.  An  estimate  for  the  maximum  gain  can  be  obtained  when  k  « 

u3  =  z[^3z  -  1  +  hr  -  1] 

v3  =  z(/3z  -  1  -  ’'kz  -  1] 

which  shows  that  z  must  be  chosen  as  large  as  possible.  Since  z  = 

It  Is  seen  that  must  be  chosen  as  small  as  possible.  For  large 

of  z,  one  obtains: 

u3  «  z/z  +  G)  =  (1 ,47z)3 
v3  *  z/z^-*^)  *  (.68z)3 
and  the  maximum  gain  will  become: 

f(u-v)  *  ^z  (1.47  -  .68)  -  .68z 

*  .34  Bq /61 

In  terms  of  A8,  It  follows  then  that: 

®max  *  A(.34  §q/ Bi *  .24djJq 


e.,  for 
*Sn* 


Bq/Bi , 
values 
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APPENDIX  II 
Eigenvalues  of  (TqSj) 

It  can  be  shown  rather  easily  that  the  4  eigenvalues  Af  of  the  transformation 
matrix  (T0SX )  have  to  fulfil  the  following  4th  order  equation: 

[A2  -  A(Zjj(1 )coshx  +  Z2 2 (1 ))  +  coshx]  •  (A2.1) 

•  [A2  -  A(Zx*(2)coshx  +  Z22(2))  +  coshx]  + 

-  slnhx2  [1  -  AZ11(1 )]  [1  -  AZi*(2)]  =  0 

with  x  =  stj 

If  one  assumes  now  that  cosh  sii  «  1 ,  one  finds  easily: 

(A  -  X,)(A  -  X2)(A  -  XjXA  -  XJ  -  s1nhx2[l  -  AZU(1)][1  -  AZ,t(2)]  «  0  (A2.2) 

where  the  X^  are  the  elgenvaleus  of  the  unperturbed  transmission  Matrix  T0> 
namely: 

Xx  ■  exp  j0,;  X2  *  exp  -j©i  (A2.3) 

X3  «  exp  J02;  X„  ■  exp  -je2 

Physically,  X,  and  X2  are  given  by  the  phase  shift  0i  a  wave  of  frequency 
uij  will  suffer  through  one  single  but  complete  section  of  the  periodic  waveguide. 
Similarly,  02  Is  the  phase  shift  a  wave  of  frequency  u>2  suffers  through  one 
single  but  complete  section  of  the  periodic  waveguide. 
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We  expect  the  eigenvalues  A-j  of  the  dispersion  equation  (A2.1)  to  be  close  to 
the  eigenvalues  Xi  of  the  unperturbed  matrix.  Therefore,  we  will  use  the 
ansatz: 


Aj  =  Xj(l  +  ej) 

The  terms  1  -  AjZxl  In  relation  (A2.2)  can  then  be  simplified: 


(A2.4) 


and  one  finds: 


1  -  AjZn  si-  XjZn 


1  -  X,Znd)  =  -^(Xj  -  X2)  (1  -  JKxjO)) 


1  -  X2Zxl(l)  = 


1  *  X ,ZX1(2) 


^■(X, 


X2)  0  +  JKxxO)) 


-  X,)  (1  +  jK„(2)) 


(A2.5) 


1  -  X4Zx*(2)  = 


~  (X,  -  X4)  (1  -  jKjX(2)) 


These  relations  show  the  Importance  of  the  value  of  Kn,  l.e. ,  Z„.  The 
real  part  of  Zxx  was  given  In  relation  (50),  namely: 


Real  Zn  *  Real  Z22  ■  cos  8 


(A2.6) 


cos  8  *  cosB*!  coskt2  -  o  slnRij  s1nkl2 
9  *  (1  +  t)/Z/c 

e  -  dielectric  constant 


Similarly,  the  Imaginary  part  of  Zn  can  be  evaluated: 


ImZ,!  ■  -ImZ22  »  -jKn  sine 

3  -  (coset j  s1nkt2  +  o  slneti  coskt2) 


Furthermore,  since  the  determinant  of  the  matrix  (Z)  Is  equal  to  1: 


Z11Z22  "  Z12Z21  a  1 


One  finds  easily  with  Zn  3  Z22  *  cose  +  Kxis1n8 


~(Kii  +  1 )  s1n2e  *  ZX2Z21  *  (£  "  ^  s1n2Btx 

4^ 


resulting  In  the  relation: 


K  a  ±j 


\2  .4.2 t 


4e  s1n2e 


e  *  dielectric  constant 


(A2.7) 


(A2.8) 


Clearly,  relations  (A2.8)  and  (A2.7)  show  that  Kxx  may  become  Infinite  as  8 
approaches  kir,  and  that  the  product 


A»  -  A, 


O  +  JK,i)  -  j  Sine  (1  +  JK,X) 


approaches  the  value  -  sine  Klx  as  6  ■+■  kir.  For  6  3  kir,  the  value  stays  generally 
finite  If  8tx  t  kir. 
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We  will  now  differentiate  between  cases. 


Case  I:  X-t  f  Xj  (1  f  j) 

In  this  case,  the  eigenvalues  of  the  unperturbed  matrix  are  all  different 
from  each  other.  Using  (A2.4)  and  (A2.5)  one  finds: 


e2  " 


es  m 


c 


-  stnh2x  [1  -  jKlt(l)] 


-  s1nh2x  [1  +  jKn(l)] 


-  s1nh2x  LI  +  jK,,(2)] 


-  s1nh2x  [1  -  jKn(2)] 


1 

-  X2Zii 

*(2) 

2(X1  - 

x,)(x, 

-  X,) 

1 

-  x2Zn 

*(2) 

2(X,  - 

x,)(x2 

-  X„) 

1 

*  XjZn 

(1) 

2(X,  - 

xjLx, 

"  x2) 

1 

*  X%Zji 

(1) 

2(X,  - 

X1)(X„ 

"  x2) 

(A2 


It  Is  seen  that  ej  Is  proportional  to  the  square  of  the  depth  of  modulation  5 
(slnhx  »  x  *  6/0! B2  Ax/2),  and  as  such,  relatively  low  gain  will  be  obtained. 
It  Is  Interesting  to  note  that  Kn  can  become  Infinite  for  sin  e  -  0.  In  this 
case,  however,  Xt  and  X2  are  nearly  equal,  so  that  the  case  of  the  it  modes 
must  be  considered. 

Case  II:  The  singular  w  mode 

Let  us  assume  that  6,  Is  nearly  equal  to  kir,  and  that  02  f  kir.  To  solve 
the  dispersion  equation,  let  us  set: 


X2  ■  X,(l  +  S2)  *  cos  8j  •  ±1 
|i2|  «  1.  f  X, 
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Let  us  solve  again  the  dispersion  equation  (A2.2)  with  the  ansatz: 


A  *  X,(l  +  et) 

The  first  term  in  equation  (A2.2)  becomes  equal  to: 

Xi2  e(e  -  S2)(Xa  -  XsKX*  -  X„) 

The  second  term  of  the  same  equation  can  be  evaluated  by  using  the  following 
relations:  (§  0) 


1  -  Xjl  +  c,)  Zn(l)  =  (1  +  El)  [1  -  XxZlx (11 )]  -  e, 


-  (1  +  ej 


- o  (Xi  ~  X2)  (1  -  jKu(l)) 


“  Ei 


=  -  XiKIJ(l)s1n01  -  e,  -  -jX,ImZn(l) 


Consequently,  with  X,  *  ±1,  one  obtains: 


e(e  -  S2) 


-i 


sinh2x 

2 


ImZnO) 


1  -  Zit(2) 

1  f  COS02 


(A2.10) 


.  -j  iljj*!*  (imZ^d)) 

For  S  0,  It  follows  that  e  Is  proportional  to  slnhx.  It  Is  Interesting  to  note, 
however,  that  the  right  hand  of  equation  (A2.7)  can  become  Infinitely  large  If 
the  denominator  1  *  cos8x  vanishes.  This  occurs  for  cosex  ■  2irk  If  0i  *  2wk' , 


±1 


l  f  cose 


r 

2 


/ 
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or  for  cose2  *  (2k  +  1  )ir  If  0j  «  (2k  ‘  +  1)ir.  Therefore,  It  Is  expected  that 
a  maximum  value  of  e  will  be  obtained  when  8X  *  02  +  2lor  and  0!  *  nir.  We 
will  consider  In  the  following  the  dual  it  modes. 

Case  III:  The  dual  it  inodes 

The  dual  ir  modes  will  be  characterized  by: 

0X  *  02  +  2kir  *  nir 

To  solve  the  equation  (A2.2),  we  will  use  the  ansatz: 

A  *  X,(l  +  et) 

X-j  -  Xj(l  +  Sj)  1  >  2 

Similar  to  the  discussion  above,  one  finds  for  the  left  term  In  equation  (A2.2): 
X/  -  S2)(et  -  S$)(ei  -  §0 

and  the  right  term  of  the  same  equation  can  be  evaluated  by  using  the  two 
relations. 


1  -  XZn(l)  =  -jA^ImZnO)) 

1  -  X(Z,*(2))  *  jXiI«,Zai(2) 

It  follows  therefore  that  the  equation  (A2.2)  becomes: 

ex(e,  -  (.Me.  -  ls)(ei  -  Ij  -  s1nh*x  WZn(l)  la(Zu(2)) 
which  results  for  perfect  synchronism  (0X  ■  0a  *  kir): 

e,2  •  ±  slnhx  /[Im  ZJ1)  1*^(2)] 


(A2.11) 


(A2.12H 


From  the  condition  that  the  determinant  of  the  matrix  (Z)  Is  equal  to  1,  namely: 


(RealfZjj))2  +  (Im(Zir))a  -  Z12Z21  *  1 

and  that  Real  (Zu)  ■  ±1  at  0  *  kir.  It  follows  that: 

Im(Zn)  =  ±  ^12Z21  =  ±  ~~  slnBtj 

2/e 

and: 

ex2  =  ±  ~~~  slnhx  [isIngjA,  s1n$2t1]/2 

Relation  (A2.13)  shows  that  ex  Is  proportional  to  the  square  root  of  slnhx  *  x 
showing  that  large  gains  can  be  obtained. 


(A2.13) 
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Equation  (39)  In  Section  5.0  required  the  synchronism  of  the  two  waves, 
I.e.,  the  wavevector  Bp  of  the  pump  wave  was  set  equal  to  the  sum  of  the  two 
wavevectors  Bx  and  $2.  Physically,  this  means  that  the  waves  are  In  complete 
synchronism  within  the  empty  waveguide  sections.  However,  under  these  con¬ 
ditions,  the  pump  wave  will  advance  during  the  passage  of  the  waves  Ft  and 
F2  through  the  dielectric  section.  Therefore,  the  second  relation  of  (39) 

Is  not  necessarily  fulfilled.  Instead,  one  must  write: 


Bp  *  B,  +  B2  +  AB 


Following  the  approval  shown  In  Section  5.0,  one  finds  then  that  the  value 
of  s  must  be  replaced  by: 


\  'f  6*6!  B2  -  AB* 


(A3.1 ) 


(A3. 2) 


1  as  derived  In  Section  4.0.  Furthermore,  it  Is  easily  shown  that  the  coupling 

j  matrix  (sx)  In  relation  (46)  must  be  replaced  by  a  matrix  (Si')  given  by: 

.  ol  0  o2  0 

*  0  10  0 

s, '  »  * 

|  as  0  Oi  0 

0  0  0  1 

[  L 

[; 
f] 

D 


with 


✓ 


a,  =  [coshs'Aj  +  j  ^.slnhs'*,]  exp  -  j  t, 

a2  =  i  |ff  slnhs'tj  exp  -  j  i,  (A4.4) 

a3  =  -j  Ma.  slnhs'i.,  exp  j  ^  t. 

2s'  2 


Consequently,  the  4th  order  equation  for  the  Eigenvalues  (see  Equation  A2.1) 
will  become: 

[A2  -  A(ZM(1)  <Xj  +  Z22(l))  +*,]  • 

[A2  -  A(Zj*(2)  a*  +  Z22(l))  +  at]  + 

"  °2a3  [1  "  AZ1X(1)]  [1  -  AZjt(2)]  =  0 


Note  that: 

a2a9  55  sinh2s%  (A3. 6) 

4S 1 

aiat  -  a2°S  =  1 

Let  us  assume  now  that  oj  *  o2*  *  1.  One  obtains  then  the  equation: 

(A  -  XJ(A  -  X2 ) (A  -  X*)(A  -  XO  «  slnh  s'l,  . 

*  Cl  -  AZ, *(!)][*  -  AZ,t(2)] 
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